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We consider a problem that models fluid flow in a 
thin domain bounded by two surfaces. One of the 
surfaces is rough and moving, whereas the other is 
flat and stationary. The problem involves two small 
parameters e and \x that describe film thickness 
and roughness wavelength, respectively Depending 
on the ratio X = e//x, three different flow regimes 
are obtained in the limit as both of them tend to 
zero. Time-dependent equations of Reynolds type 
are obtained in all three cases (Stokes roughness, 
Reynolds roughness and high-frequency roughness 
regime). The derivations of the limiting equations are 
based on formal expansions in the parameters e and \i. 



1. Introduction 

The fundamental problem in lubrication theory is to 
describe fluid flow in a gap between two adjacent 
surfaces which are in relative motion. In the 
incompressible case, the main unknown is the pressure 
of the fluid. Having resolved the pressure it is possible 
to compute other fundamental quantities such as 
the velocity field and the forces on the bounding 
surfaces. To increase the hydrodynamic performance 
in various lubricated machine elements, for example 
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journal bearings and thrust bearings, it is important to understand the influence of surface 
roughness. In this connection, one encounters various approaches, commonly based on the 
equation proposed by Osborne Reynolds in 1886 [1]. Although a number of averaging 
methods considering surfaces roughness have been proposed over the last 40 years (e.g. [2-4]), 
homogenization has prevailed as the proper way to average [5,6]. Homogenization is a rigorous 
mathematical theory that takes into account information about local effects on the microscopic 
level [7]. 

This study is concerned with the asymptotic behaviour of Stokes flow in a narrow gap 
described by two small parameters e and \i. The parameter e is related to the distance between the 
surface, whereas \i is the wavelength of the periodic roughness. In many problems involving two 
small parameters, the way in which the parameters tend to zero is primordial and the limiting 
equations may be different whether e tends to zero faster, slower or at the same rate as \x. Using 
formal asymptotic expansions in the evolution Stokes equations, we show that three different 
asymptotic solutions, i.e. three different flow regimes, exist in the limit as e > 0 and \i > 0 tend to 
zero depending on whether the limiting ratio 

k = lim — 

(6,/x)^(0,0) [I 

equals zero, a positive number or oo. In all three flow regimes, the limiting pressure is governed 
by a two-dimensional equation of Reynolds type whose coefficients take into account the fine 
microstructure of the surface, i.e. a homogenized equation. The situation can be summarized 
as follows: 



Stokes roughness regime. The case when 0 < k < oo. One finds that the coefficients of the 
homogenized equation are obtained by solving three-dimensional so-called cell problems 
which depend on the parameter k. 

Reynolds roughness regime. The case when k = 0. The cell problems are two-dimensional 
and the proposed averaged equation appears in, for example, [3,8,9]. The same limiting 
equations are obtained if one lets k -> 0 in the Stokes roughness. 

High-frequency roughness regime. The case when k = oo. We obtain a limiting equation of 
very easy and cheap treatment. The same limiting equations are obtained if one lets k -> 
oo in the Stokes roughness. 

This work is closely related to the studies by Bayada & Chambat [4,10] and Benhaboucha et al. 
[11], who considered the stationary case, i.e. only the flat surface is moving. The main novelty 
is the treatment of the unstationary case (the rough surface is moving) as well as the way that 
e and \x tend to zero. The paper is organized as follows: §2 is devoted to the formulation of the 
problem and basic notations. Section 3 contains a summary of the main results of this work. In 
§4, we define the formal asymptotic expansions, the corresponding change of variables, domains 
and differential operators for the problem. Section 5 is concerned with the Stokes roughness, with 
constant ratio €/fi = k. This is the case analysed in [4,10]. Section 6 is devoted to the case e = /x 2 , 
which corresponds to Reynolds roughness. We apply the asymptotic expansion method in one 
parameter and derive the homogenized Reynolds equations. The last section deals with the case 
fx = e 2 , which belongs to the high-frequency roughness regime. We obtain the classical Reynolds 
equations with truncated film thickness. We note that neither e = fi 2 nor fi = e 2 is covered in 
[4,10], whereas [11] only covers e = n 2 . Evidently, as mentioned in [4,10], identical equations 
are obtained if one lets k -> 0 and k — > oo in the Stokes roughness regime. However, from a 
mathematical point of view, there is no apparent reason why taking limits in such different ways 
would yield the same result. 

For clarity, the main results are presented as 'theorems' and their derivations as 'proofs', 
although the method of formal expansion is not rigorous by mathematical standards. Choosing 
this style, we hope to make the paper accessible to a wider audience. We stress, however, that all 
calculations (including limit processes) can be made rigorous. 



2. Problem formulation and basic notations 



This study is concerned with thin film hydrodynamic lubrication of rough surfaces. For simplicity, 
we suppose that one of the surfaces is rough and moves with velocity v = (v\, v 2 , 0) and that the 
other is flat and stationary As the rough surface is moving, the film thickness varies in both 
space and time, thus rendering the problem unstationary A point in space (R 3 ) is denoted as 
x = (xi,x 2 ,x 3 ) / and fisa time variable that belongs to the interval [0, T]. The problem considered 
is the evolution Stokes system 

du 

vAw + Vp = 0 (2.1) 



dt 



and 



div u = 0, 



(2.2) 



where v (viscosity) is a constant, and u — (u\,u 2r u 3 ) (velocity field) and p (pressure) are unknown. 
We shall write 



x f = (x 1 ,x 2 ), 



X3 



y 



and r : 



t 



(2.3) 



where e and \x are two small parameters. The basic idea of the homogenization method is to treat 
x' ' ,y' ' , t and x as independent variables. Equations (2.1) and (2.2) are assumed to hold in a moving 
space domain fi €lJL {t), defined by 

Q €lJi {t)= jx = (xi,x 2 ,x 3 )£R 3 , x f = (x 1 ,x 2 )ecD, 0<x 3 < eH (x f , ^-,t, -0 J, 

where co is an open connected set in R 2 with smooth boundary, outward unit normally denoted 
by h and the function H{x' , y' ', t, r) describes the geometry of the upper surface. H is assumed to be 
Y-periodic in y' , Y = [0, 1] x [0, 1] being the cell of periodicity and T-periodic in r. More precisely, 

H{x f ,y f , t, x) = h 0 (x f - tv) + h per (y f - xv), 

where /zo describes the global film thickness, whereas the Y-periodic function /z per represents the 
roughness. Thus, e is related to the film thickness, whereas \i is the wavelength of the roughness. 
Moreover, we define the 'minimum film thickness' 

HV, t) = h 0 (x f - tv) + miny G y/z p er(y / )- 
The boundaries of £2 €lJb {t) are 

U^(t)=^xeR 3 :x f = (x 1/ x 2 )ea) / x 3 = €H (^x f ,^,t, ^ •, 

27 € - (*) = {x e R 3 : x' = {x lr x 2 ) e co, x 3 = 0} 
and E^{t) = |x e R 3 : x' = (x lf x 2 ) e dco, 0<x 3 <eH (x f , t, J . 

We assume the following no-slip boundary conditions: 



dH 

V1,V2,€\—+V' VH 



u = 0 



and 



on (t), 
on S-^t) 



where g = (gi,gi,g3) is some given function and the initial condition 



el4(xi,x 2 , X j) 



where Vff = (8H/dx 1 , dH/dx 2 , 0). 

For convenience, we use the notation 



on^ f/1 (0) x (0), 



fdy', r 



fdz and f z 



fdz 



for integrals of a function/. Moreover, we denote by e\, e 2 , £3 the standard basis vectors in R 3 . 

Finally, to ensure the existence of u, we must require some compatibility between the boundary 
conditions and H. To this end, it is assumed that g is a C 1 vector field defined on R 3 such that 
divg = 0, g(x 1 ,x 2 , 0) = (0, 0, 0), 

dH 

gi(x 1/ x 2 ,H) = vi, g2(xi,x 2 ,H) = v 2 and g3(x 1 ,x 2 ,H) = — - + v • VH 

ot 

and 



(2.4) 



r dH f 




— dx f + 




L dt J 


dco V 



gdz\ • n dS = 0 



for all (y',£,r). 



3. Formal asymptotic expansion in e and /x 

We analyse the asymptotic behaviour of the equations of motion (2.1) and (2.2). We define the 
following expansions for u and p: 

oo oo 

u ( Xft ) = J2 e n n m u n > m (x',z,y',t,r) 

n=0 m=0 

and 

oo oo 

p(x,t)= ^VV'^UyUr), 

n=-2m=0 

where x' , z, y' and x are defined by (2.3) though subsequently treated as independent variables. As 
the roughness is periodic, it is assumed that u n,m {x f ,z,y f , t, x) and p n,m {x f ,z,y f , t, x) are Y-periodic 
in y' and T-periodic in r. It is convenient to define also the following domains: 



V{y f r t r x) 


= {(x',z)ecox R 


0<z<H{x',y f , t,x)}, 


%{x',t,x) 


= {{y',z)eY x R 


0<z<H(x f ,y',t,x)}, 




= {(y f ,z)eYx R 


H*(xf ',t)<z<H(x' \y' ',t,x)} r 


Q*(t) 


= {(x',z)eco x R 


0<z<H*(x',t)}, 


Z*(x',t) 


= {(y f ,z)eYx R 


0<z<H*(xf,t)} 


Y z {x f ,t,x) 


= {y f eY:(y , ,z)e^(x / ,t,x)}. 



The boundaries of Q(y' ' , t, r) and Q*(t) are 

E + {y r , t, r) = {{x',z) ecoxR:z = H(x' ' ,y' ', t, r)}, 
E ~ {y' r t,r) = {{x', z) e oj x R : z = 0}, 
^(y', f, r) = {(x',z) edcuxR:0<z< H(x' ' ,y' ', t, r)}, 
= {(x', z)ewxlR:z = H*(x', £)}, 
17*-(0 = {(^ / z)e(»xR:z = 0} 
and 27*^(0 = {{x f , z) e dco x R : 0 < z < H*(x', f)}. 

Note that B(x', f, r) and B*(x', £) do not have lateral boundaries because of the periodicity of Y. The 
boundary of Y z is denoted by 

dY z {x f , t,x) = {y f eY:z = H(x' ' ,y' ', t, r)}. 

(a) Differential operators 

Assume u(x, t) = v(x! ,z,y' , t, r). Then, 

du dv 1 dv 
dt dt /x dr ' 



( dv dv \ 1 dv dv \ l / n 

Vu = I , ,0 + - , ,0 + - 0,0, 

\dx! dx 2 J [i \dy\ dy 2 ) e\ 



. d 2 d 2 \ 2 d 2 d 2 \ i d 2 d 2 . 

and Au = — = H ~ v H 1 iH ~ — ~ H ~ u + ^=r — n v. 

dx 2 dx\) ii\dx x dy x dx 2 dy 2 ) /x 2 \dy 2 dy 2 ) e 2 \dz 2 * 



Define 



and 



Moreover, 



d d \ / a a \ t ( d 

V* = [—,—,0), Vy = __0 and V z = 0,0,- 

d 2 d 2 d 2 d 2 , a 2 

= — A v / = — y and A z = — r-. 

ax 2 ax 2 y ay 2 ay 2 az 2 



div x /i; = V X ' - v and dhyv = • v. 



4. Stokes roughness 



Figure 1 describes the case when € and /x tend to zero with constant ratio 0 < k < oo. That is, 
we assume that fi = €/k, where A is a positive constant. We define the following asymptotic 
expansions: 

oo 

u(x, t) = J2 £ n u n {x', z, y', t, r) (4.1) 

n=0 

and 

oo 

p(x,t) = J2 *W,z,y',t,r)- (4-2) 

n=-2 

Inserting (4.1) and (4.2) into (2.1) and equating terms of the same order using 

a a a a 

dt ~* dt + € dr' 

X 1 

V* + -Vy + -V z 



global film thickness 

■ rough and moving surface 

■ flat and stationary surface 




Figure 1. Stokes roughness (/x = e /k). (Online version in colour.) 



gives 



: kVy>p + V z p~ 2 = 0, 



-=■ : -vk 2 Ay u° - vA z u° + V x >p~ 2 + kS/y,p~ x + V z p~ l = 0 



and 



1 3m° 



2 7/ 0 



3^ 



2 7/ 0 



:A — y I 2A— ^ \-2k—^ h A z w 2 + X 2 A V 'U 1 

e dx \ dxxdyi dx 2 dy 2 y 



+ Vx'p- 1 + XVy'p 0 + V Z ?° = 0. 

Similarly for (2.2), we have 



and 



1 n aw§ 

- :Adhyw° + =0 

6 * dz 



du 

e° : div x /w° + k divyu 1 + — ^ = 0. 

y 3z 



The boundary conditions are 



o 



and 



(gl/g2,0) 

0,0, — + u • V x /H ) 



on i: - 
on £ w 



0 

(0,0, g 3 ) 



on ir 
on 



and initial conditions are 



and 



6° : u° = (LI?, 17°, 0) on Q x {f = 0} 
e 1 =(0,0, U%) on^2 x {f = 0}. 



(a) Analysis of equations 

The main result pertaining to the Stokes roughness is as follows. 
Theorem 4.1. The leading term u° in expansion (4.1) for u is given by 

i=\ dXi 

where a l (i = 0, 1, 2) is a solution of the periodic cell problems (4.12) and (4.13) and £/ze leading term p~ 2 in 
expansion (4.2) for e 2 p is a solution of the boundary value problem 

div x /(A A V x /j^- 2 + &*) + =0 in cox (0, T] (4.9) 

and 

{A k V x ,p- 2 + b x -g* y ')-h = 0 onda)x(0,T], (4.10) 

where A k and b x are calculated by (4.17) and (4.18), respectively. 
Proof If we write (4.3) in component form, we obtain 

=0, J —=0 and ^— =0, 
3yi 3i/2 9z 

hence p~ 2 {x f ,y f ,z,t, x) =p _2 (x / ,f,r). We are looking for solutions w° of the form (4.8) and p~ x of 
the form 

where a 1 = a l (x! ,z,y' ,t, x) and q l = q l (x' ,z,y r ,t, x) are to be determined. Clearly, (4.8) and (4.11) 
satisfy (4.4) and (4.6) if 

v^Aya* + vA z a j = Wyq* + V z if + e { in B (z = 0, 1,2) (4.12) 

and 

k divya 1 + = 0 in B (z = 0, 1, 2). (4.13) 

J dz 

The above systems of equations are called cell problems, whose solutions a 1 and q l are Y-periodic, 
and the boundary conditions are 

(v^^+v^h)), (/ = 0) onS+ 

0, (f = l,2) 

and a z ' = 0, 0 = 0,1,2) onS". 

It can be verified that each a 1 is uniquely determined by (4.12) and (4.13). 



Multiplying (4.7) by (j)(x r ) e C 1 (<w) and integrating by parts using the Gauss-Green theorem, 
we obtain 



0 = 



n 1 du \\ 

div x /w u + k div y >u l + J <p{x f ) dz dy' dx' 



Y J 
+ 



y J 



y J 
+ 



divj'i/°0 dx f dz dy' + 

-u° • V J /0dzdx / dy' + 
0g • ndS(x' ,z) dy + 
-u° Z ■ V x /0dy' dx f + 



X divy'U 1 H I 0 dy' dz dx' 



3z 



y J 



9B 



0(^,w 2 , ,^)-ndS(x / ,z)dy / 



0A(t/{,^,O)-ftdS(y / ,z)dx / 



3co - 



0g • (hi, fi2, 0) dS^, z) dy 



-v • V x /H0 + + u • V x /f/) 0 dx' dy 
0g zi/ -ndS(^). 

As 0 is arbitrary, it holds that 



y J 



div x >u° + =0 in a; x (0, T] 

3f 



and 



u° -g^)-h = 0 ondcox(0,T]. 



By integrating (4.8), we obtain 



where 



A 



a| 0/ 



and 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



Inserting (4.16) into (4.14) and (4.15), we obtain the homogenized Reynolds equation (4.9) with 
the boundary (4.10). ■ 



5. Reynolds roughness 



9 



Figure 2 describes the case when the wavelength of the roughness is much greater than the film 
thickness, i.e. /x ^ 6. This case can be studied by assuming that e is a function of /x such that 

lim^=0. 

/l^O /X 



For simplicity, we shall assume that e = /x 2 . We postulate the following expansions for u and p : 



u(x, t) = J2 ^ n u n {x' f z,y', t, r) (5.1) 

n=0 



and 



p(x,t)= H n ftf,z,l/,t,r). (5.2) 

n=-4 



Plugging (5.1) and (5.2) into (2.1) and equating terms of the same order using 



a a l a 

dt dt n dr ' 

1 1 
/x /x 2 

2 / a 2 a 2 \ i l 

and A -> A x / H I — — — + — — — + -^Ay + — A z 

/X \ dx\dy\ dx 2 dy 2 l /x /x 



gives 



1 

"76 



/x' 
1 

75 



V z p~ 4 = 0, (5.3) 
Vyjy- 4 + V z p- 3 = 0 (5.4) 



/x 

i a 2 i/° 



and : -v^- + V *^~ 4 + V y^" 3 + V ^" 2 = 0- (5-5) 

/x* dz z 

Similarly for (2.2), we have 

i aw° 

^:-^=0, (5.6) 

1 n 3wJ 

:divy>u° + —^=0 (5.7) 



/x y az 



dul 



and /x°:div x /W° + div 1/ /W 1 + ^=0. (5.8) 

* az 



Figure 2. Reynolds roughness (jjl ^> e). (Online version in colour.) 



The boundaries conditions are 



H°:u°: 



(vi, v 2 ,0) onE+ 
0 on E~ 

dH 

0,0, — + u-V l/ H) on 17 



on E~ 
on E w 



and 



and initial conditions are 



2 2 



0,0, + V x /H) 



0 

(0,0, £ 3 ) 



on E~ 
on £ w 



fi°:u° = (u\, u\, 0) on (J ^ /0//x x {j/} x {0}, 
and /x 2 :w 2 = (0,0,^). 



b 



?3 



(a) Analysis of equations 

The main result is as follows. 

Theorem 5.1. The leading term u° in expansion (5.1) for u is given by 

n z(z - H) . _ 4 _ _« z 

2v ( x ' v +v y' v ) + H V ' 

where p~ 4 , the leading term in expansion {5.2) for \i^p, is a solution of the boundary value problem 



div x /(-AV x ^- 4 + b) + — — =0 in co 

dt 



and 



(-AV x >p~ 4 + b-g zV )-h = 0 ondco, 
where A and b are given by (5.22) and (5.23), and 



-3 ^ 3P" 4 
i=l M ' 



(5.9) 

(5.10) 
(5.11) 

(5.12) 



where qi (i = 0, 1, 2) is a periodic solution of the cell problems (5.16) and (5.17). 



Proof. Equations (5.3) and (5.4) say that 



3r?- 4 n 9r?- 4 n dp~ A n , a?" 3 „ 
r -0, =0, =0 and =0. 



az 



3yi 



dy 2 



dz 



Therefore, 



p~ 4 (x f , z,y', t,r) = p~^{x', t, r) and p~ 3 (x' ' ,z,y' ', t,r) = p~ 3 {x! ,y' , t,r). 



From (5.6), we deduce u® = 0 in Q because of the boundary conditions for u®. Thus, (5.5) in 
component form becomes 

d 2 u\ dv~ A dp~ 3 
az 2 axi 3yi 
d 2 u° 7 dv~ 4 dp~ 3 
az 2 ax 2 ay 2 

9 r 2 _ n 



and 



az 



Hence, the first two equations may be written as 

d 2 u 0 



dz 2 



+ v x ^- 4 + v^- 3 = o. 



(5.13) 



Integrating (5.13) with respect to z and taking into account the boundary values of u°, we get (5.9). 
Integrating (5.9) once more, we obtain 



Multiplying (5.7) with 0(y') and integrating over B using the Gauss-Green theorem gives 

dul 



(5.14) 



div^ o + ^j0(yOdy'dz 
-u° • V y /0dy / dz+ (j){u\,u\ f u\) •ndS(y / ,z) 



y J 
+ 



V y /0dzdy' 



0 vi 



an dH 

dyi dy 2 ' 



-m° •V / 0+—0dy / 



( div *' M ° 2+ f) 



Udy', 



for all smooth and Y-periodic 0. Hence, 



,z . an 



divyu 0 + — =0, in Y. 
y dz 



(5.15) 



Inserting (5.14) into (5.15), it is seen that p 3 can be written in the form (5.12), where qi is periodic 
solutions of 

H 3 



div y' " 7T" ( y y"?i + e i) = °' in Y (' = !' 2 ) 



and 



ff 3 



(5.16) 



H \ 3H 



This implies 



u u = 

12v 



dyi dyi 

dqi 1 8^2 
V 9y 2 9y 2 / 



H 



Multiplying (5.8) with (p(x f ) e C 1 (&») and integrating gives 
0 = 



y J 



/ n 1 

div x /w° + dhy w 1 + — ^ 0(x') dz dy dx r 
o V 9z ' 



y J 
+ 

y J 
+ 



div x /w°0dy / dx' + 



-u° • V x /0dzdx' dy' + 



dui 



divx'U 1 + ) 0 dx' dy' dz 



<t>(u\,u\,vi§ • hdS{x' r z) dy' 



-u" • V^dzdx'dy' + 



y J 



(pg • hdS(x f ,z) dy' 



v\, V2, h v • V x /H ) • n dS(x',z) dy' 



-u° •V x /0dx / dy / + 



y J 



(-v • V x 



3oj J 



0g • hdS(x f ,z)dy r 

-U° • V x /0+ — — 0dx' + 



-ndS^z) 



As 0 is arbitrary, it holds that 



and 



where 



div x /W u H =0 m 60 

dt 



=^ y —y' 

(u° -g zy )-n = 0 ondco, 



— prZ" A 



and 



y 12v 





dyi 




1 + — 

8y 2 


V o 


0 





dy' 



(5.18) 



o 



b 



5 



(5.19) 
(5.20) 
(5.21) 

(5.22) 



(5.23) 



Inserting (5.21) into (5.19) and (5.20), we obtain the homogenized Reynolds equation (5.10) with 
boundary condition (5.11). ■ 

Remark 5.2. The equations of theorem 5.1 can also be obtained by letting X —*■ 0 in the equations 
of theorem 4.1. 



6. High-frequency roughness regime 



Figure 3 illustrates the case when the roughness wavelength is small compared with the film 
thickness, i.e. e ^> \i. To this end, one assumes that \i is a function of e such that 

lim — ^— = oo. 

For simplicity, we shall assume /x = e 2 . We postulate the following expansions: 

oo 

u(x, t) = J2 £ n u n {x', z, y', t, r) (6.1) 

and 



n=0 



p(x,t)= e n V n {x r ,z,y\t,T). 

n=-2 

Plug (6.1) and (6.2) into (2.1) and equating the terms of the same order using 

a a 1 a 

i l 

v^v x > + ^v / + -v z 



and 

that gives 



2 / a 2 a 2 \ 1 l 



(6.2) 



and 



? :-vAj/ M 0 + V y p =0, 
T : -vAyU 1 + V z p~ 2 + Vy^ -1 = 0 



e 3- 

1 3m° 



v 2 



2„0 



a 2 i/' 



+ 2- 



a 2 u° aV 



6 2 ' dx " y dX\dy\ ' " dX2dy2 ' 3z 2 

+ v*p- 2 + v z p- 1 + v 3 ,p° = o. 

Similarly for (2.2), we have 



1 



■ :div 1 /«° = 0 



i i 

- idivyu 1 + — A = 0 



and 



e° : div X 'W u + div V 'U z + 



,,0 



a^ 
az 



= o. 



The boundary conditions are 



and 



0 on E~ 

dH 

0,0, — + v-V x >H) onS+ 



0 

(0,0,g 3 ) 



on E~ 
on r 1 " 



(6.3) 
(6.4) 

(6.5) 

(6.6) 
(6.7) 
(6.8) 



global film thickness 

- rough and moving surface 

- flat and stationary surface 




Figure 3. High-frequency roughness regime (e ^ fi). (Online version in colour.) 



and with initial conditions 



and 



6°:i/° = (U?,U°,0) 
e l :u l = (0,0,17°). 



(a) Analysis of equations 



The main result is as follows. 

Theorem 6.1. The leading term u° in expansion (6.1) for u is given by 



z(z - H* 
2v 



(6.9) 



where p 2 , the leading term in expansion (6.2) for e 2 p, is a solution of the boundary value problem 



H 



*3 



div x / ( — : Vv'P z + : — v 1 + — — =0 in co 

12v xy 2 dt 



3H* 



and 



H 



*3 



H* 



, v ^x'P 2 + t ^-v-g z *)-n = 0 ondco. 



(6.10) 



(6.11) 



Proof. Multiplying (6.3) with 0(i/) = (01/02/03) such that 4>{y f ) vanishes near 9Y Z and 
integrating over Y z gives 



0: 



-v Ay>u° ■ 0 + V yf p~ 2 ■ (p dy r 



v ^ V y /W ° • v y'0z ~V 2 div y /0 dy'. 

: i=l 



If z < H*, then Y z = Y and 3 Y z = 0 so we can choose 0 = u°. By (6.6), we obtain 



0 = 



Y^\Vy>u?\ 2 dy f ^vVy>u 0 = 0 inY; 



z=l 



if z > H*, then we can choose 0 Z - = — Uj. Thus, 



0 = 



^IV^pdy' = 0=^tt? = «i in Y z (1 = 1,2). 



In both cases, Vy/jt? 2 = 0. Summing up, it holds that 

Vyu? = 0 inY ifz<H* 0 = 1,2,3), 
u\ = v { inY z if z>H* (z = l,2) 



(6.12a) 
(6.12b) 



and 



V y >p~ z = 0 in I 



Integrating (6.7) and using the Gauss-Green theorem gives 

i , 

divy/w 1 H dy dz 



dz 



u\,u\,u^) • ndS 



9B+ 



K,^,^)| z=z ..(0,0,-l)dy / 



= -^(x / ,z / ,f,r)|Y z /|. 
Thus, w° = 0 for all 0 < z' < H. Thus, (6.4) and 6.7 become 



y Ay w 1 + Vyp- 1 + V z p~ 2 = 0 in I 



and 



divy'U 1 =0 in I 



(6.12c) 



(6.13a) 
(6.132?) 



with the boundary conditions 

/ da \ 

0,0, — + u-V x /H) onS+ 
0 onS~. 
Multiplying (6.13a) with </>(y f ' ,z) = (0i, 02, 03) vanishing on S + U S - and integrating over B gives 
0 = j (-v Ayu 1 + Vyp- 1 + V z p~ 2 ^j • 0 dy' dz 

3 

v XI V W " v y'^' dh y^ p -2 -r— d y' dz - 



z=l 



Choosing 0 = (uj, w^, 0) gives 



b 



0 = 



v^2\V y >u}\ 2 dy'dz^ Vy>u] =0 in B, (f = 1,2). 



z'=l 



Inserting this into (6.13fl) yields Vyp 1 = 0 in B. The third component of (6.13a) reads 

dp~ 



vA y >ul + — — =0 inB 
y 6 dz 



and 



dH 

~dt 
0 



+ v • V x /H on 



on X" 



Assume z < H*. Integrating (6.14) over Y gives 
0 = 



y A^ul + ^— dy' = ^— |Y| = 

y 3 9z y ^ 



3z 



~9z~ 



= 0 ifO <z<H*. 



Hence, 



-vAy>u\ = 0 inY if0<z<H*. 



From this we deduce Vyw^ = 0. Summing up if 0 < z < H* gives 

dp~ 2 



dz 



= 0, if 0 < z < H* 



and 



V y /i4 = 0, if 0<z<H* 



(6.14) 



Assume z > H*. As u° = (v\, vi, 0) in H* < z < H, we deduce from (6.8) that 



du 



divy'U 2 H -=0 inB+*. 



So if z'>H*, then 



0 = 



dz 

9 du\ , 

div y >u A H dy dz 

dz 

(u\, u\, u\) • h dS(y' , z) 

3H 3H 



(u\, u\, u\) 



dyi dy 2 ' 



l) <V 



K4 w 3)lz=z-(0,0,-l)dy / 



3H 



+ v 'V x ,H-u\\ z=z >&y f . 



Hence, 



— + v- V x ,H)dy f = 0 for all z>H*. 



Define w = u^ — (8H/dt + v • V x /H). Then, w satisfies 

dp~ 2 



-vAyW + 

J dz 



= 0 in Y z 



and w = 0 in 9Y Z 

for all z > H*. Multiplying (6.16) with w and integrating by parts gives 



0 = 



v\W y tw\ 2 + — — wdy' 

Y z dz 

v\V y >w\ 2 &y f 



asp 2 does not depend on y' and 



Y z 



w dy' = 0 



by (6.15). Consequently, Vy>w = 0 in B^, which implies 



8H 
~dt 



+ v-V x >H inBt*. 



(6.15) 



(6.16) 



Assume from now that 0 < z < H*, equations (6.12a) and (6.12fr) imply that u° does not depend 
on y' , and (6.5) reduces to 



du° I d 2 U° ?\ 9 1 n 

— - v I + A y w 2 I + V x ,p- 2 + Vz/7" 1 + V y p° = 0. 

Integrating (6.17) over Y using also that neither p~ 2 nor p _1 depends on y' gives 
du° d 2 U° 

-5--V— x +V x ^- 2 + V z ^- 1 =0 if0<z<H*. 
dr dz z 



(6.17) 



As Mg = 0, this is equivalent to 



^77° Fl 2 1J° 

T-- V TT+^P~ 2 = ° ^0<z<H* 



(6.18) 



b 



and 



dp- 1 
dz 



= 0 ifO<z<H*. 



(6.19) 



Similarly, integrating (6.8) over Y gives 
0 = 

Hence, 



n o 9wl , / n du\ 

div x >u° + dhy w 2 + - 1 dy' = I Y| div^u 0 + - 1 
y * 9z \ 3z 



div x /w°H -=0 inJ2*. 

dz 



(6.20) 



Consider now (6.18) and (6.20) as a system for (u°, u\,p 2 ). The stationary solution of this system 
is denoted by (tfi, u\,p~ 2 ), which satisfies 



b 
(3 



+ V x /p _2 = 0 in^2* 



and 



div^iTH -=0 inJ2* 

dz 



together with the boundary conditions 



dH* 



0 

(gl,g2,g3) 



oni:*" 
on E* w . 



We want to show that (u Q ,u\,p 2 ) = (w°, u\, p 2 ). To this end, define 



and 

which satisfy 
and 



u\ — u \- u\ 

~-2 -2 --2 
P =p —p , 



du 
~dx 



2Z0 



d L x£ 



+ V x /p _2 = 0 in^2* 



n 9l 4 
div x /W u H -=0 inJ2*, 

dz 



5 



(6.21) 



(6.22) 



where ir , and p are periodic in r and satisfy the boundary conditions 



(u\, u\, u L 3 ) = (0, 0, 0) on U^"U H* w . 

Multiplying (6.21) with 0 = (01,02/03) vanishing on U U and integrating over Q* 
gives 



0: 



8w° d 2 u° 



z-2 



> dx' dz 



' difi_ 
q* dr 



du° 90 



— v • p div x /0 dx r dz. 

dz dz v Y 



Choosing 0 = (wj, ^3) and using (6.22) gives 



0 = 



aw 1 

2* 9r 



+ p- 2 ^dx' dz. 
3z 



Asp 2 does not depend on z, 



^^dx'dZ: 



P- 2 



Consequently, 



1 a 

2 3r 



\u°\ 2 dx f dz = - 



dz 



^dzdx' = Q. 



dz 



dx f dz. 



Integrating this equality from r = 0 to T using the periodicity gives 

2 



0 = 











V 




.0 i 




dz 



dx f dz =0. 



Hence, w° = 0, which means u° = u°. It follows that p 2 = p 2 and = fig, so w° and p 2 are 
independent of r. Hence, equation (6.17) reduces to 



2 1# 0 



' 9z 2 



+ V x ,p~ z = 0 inJ2*, 



(6.23) 



with the boundary conditions 



(vi,v 2 ,0) onr + 
0 on Z*~ 

(gl,gl,g3) on^. 

Integrating with respect to z gives (6.9). Integrating again (6.9) with respect to z gives 



12v xp 2 



In similar fashion as above (6.20) gives 



— z* . 9H* 



div x /W° + ^-=0 in&> 



and 



(u° Z * -g z *).h = 0 on 



(6.24) 



(6.25) 



(6.26) 



Inserting (6.24) into (6.25) and (6.26), we obtain the classical Reynolds equation (6.10), in 
'minimum' film thickness H*, with the boundary condition (6.11). ■ 

Remark 6.2. The equations of theorem 6.1 can also be obtained by letting X -> oo in the 
equations of theorem 4.1. 



7. Summary and conclusion 

As mentioned in the Introduction, three flow regimes have been identified, 

(1) Stokes roughness, 0 < X < oo (see theorem 4.1), 

(2) Reynolds roughness, X = 0 (see theorem 5.1), and 

(3) high-frequency roughness regime, k = oo (see theorem 6.1). 



In all three cases, we end up with a two-dimensional equation of the form 



&iv{A k Vp + b x 
{A x Vp + b x - 



) + 



dho 
8t 



0 in a; x (0, T] 



and 




0 on d(ox(p,T], 



where p is the unknown pressure, V = (8/8xi, 8/8x2) and 



div = V-, 




«21 a 22 



an an 



and b x = 




The matrix A 1 and vector b l are macroscopic quantities known as 'flow factors'. They are 
calculated by solving local problems on a periodic cell, thus taking into account the local 
geometry, i.e. the roughness, of the problem. The expression A x Vp + b x comes from averaging 

the first two components of the velocity field u 2 ^ . As the flow is governed by an equation which 
is a generalized form of the Reynolds equation, one can say that the thin film approximation is 
valid on the macroscopic scale in all three cases. In the Stokes roughness regime, the thin film 
approximation is not valid on the microscopic scale — the local problems are periodic analogues 
of the Stokes equation and three dimensional. Consequently, the calculation of the flow factors 
comes at a high cost. However, as X tends to zero, the solution of the local problems asymptotically 
satisfies problems that are local variants of the classical Reynolds equation. Thus, the thin film 
approximation is valid also on the microscopic level in the Reynolds roughness regime. In 
conclusion, one can say that some information about the flow on the microscopic level is lost at 
the extreme cases k = 0 and X = 00. In fact, A 00 and b°° retain no information about the roughness 
(except the minimum height) and the cell problems have trivial solutions. The limiting equation 
in the high-frequency regime is exactly the classical Reynolds equation, which has been well 
studied. It can be interpreted as though the flow is prevented from entering the thin valleys of 
the rough surfaces. The information loss in the case X = 0 is due to the thin film approximation 
on the microscopic level. It would be interesting to compare A k , b x to A 0 , b° for small values of 
A as well as the corresponding flow patterns. We hope to accomplish such a study in the future 
including numerical simulations. 

As to previous studies, the present result reduces to the stationary case when dh^/dt = 0 and 
8hp er /8r = 0. Compare with eqn (17) in [4] and theorem 3.1 in [11] for the Stokes roughness; eqn 
(25) in [4] and theorem 3.2 in [11] for the Reynolds roughness; and eqn (20) in [4] for the high- 
frequency roughness. Note that, in the unstationary case, time plays only the role of a parameter 
in all three limiting equations. Although the original equation (2.1) contains the term du/dt, the 
time derivative of the unknown solution does not appear in the limiting equation nor in the 
local problems. 
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